Surveying the Nucleon-Nucleon Momentum Correlation Function in the Framework 

of Quantum Molecular Dynamics Model 
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Momentum correlation functions of the nucleon-nucleon pairs are presented for reactions with C 
isotopes bombarding a 12 C target within the framework of the isospin-dependent quantum molecular 
dynamics model. The binding-energy dependence of the momentum correlation functions is also ex- 
plored, and other factors that have an influence on momentum correlation functions are investigated. 
These factors include momentum-dependent nuclear equation of state, in-medium nucleon-nucleon 
cross sections, impact parameters, total pair momenta, and beam energy. In particular, the rise and 
the fall of the strength of momentum correlation functions at lower relative momentum are shown 
with an increase in beam energy. 

PACS numbers: 25.10.+S, 25.70.Mn, 21.45.+V, 27.20.+n 



I. INTRODUCTION 

The method of two-particle intensity interferometry 
was developed by Hanbury-Brown and Twiss (lj in the 
early 1950s. Originally they applied two photon cor- 
relation to measure the angular diameter of stars and 
other astronomical objects. Initially, the method did not 
receive universal acceptance before a number of terres- 
trial experiments were performed to confirm it. Now the 
method of intensity interferometry is commonly referred 
as the Hanbury-Brown/Twiss (HBT) effect. Although 
the original application of the HBT effect used photons 
as the detected particles, it was rapidly realized that the 
approach can be generalized to include correlation mea- 
surements for other bosons and fermions as well. The first 
measurements of the HBT effect in subatomic physics 
came from elementary-particle reactions. Goldhaber et 
al. extracted the spatial extent of the annihilation fireball 
in proton-antiproton reactions from two-pion correlations 
0- In fact, the method explores that identical particles 
situated nearby in phase-space experience quantum sta- 
tistical effects resulting from the (anti)symmetrization of 
the multi-particle wave function. For Bosons, therefore, 
the two-particle coincidence rate shows an enhancement 
at small momentum difference between the particles. The 
momentum range of this enhancement can be related to 
the size of the particle source in coordinate space. Re- 
cently, there are not only substantial experimental lit- 
eratures on the technical applications, but also a large 
number of theoretical papers on momentum correlations 
with different models from low energy to relativistic en- 
ergy heavy ion collisions. For reviews, see the references 
Ull H El More recently, HBT method has been ex- 
tended to others fields, for instance, the analogous corre- 
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lations in semiconductors and in free space aiming at the 
fermionic statistics of electrons 0, Q • 

Many experimental measurements of HBT effect have 
been performed for heavy ion collisions (HIC) at interme- 
diate energy in recent years. With the increasing of beam 
energies, nucleon-nucleon collision plays a dominate role 
during the reaction process in the intermediate energy, 
which results in an increasing importance of earlier par- 
ticle emission. By the application of the two-particle cor- 
relation function, one could obtain the information on 
particle emission and collision dynamics. 

In this energy domain, most correlation measurements 
focus on two-proton correlation functions. The shape 
of a two-proton correlation function reflects the com- 
bined effects of the Pauli blocking principle, Coulomb 
interaction and proton-proton (p — p) nuclear interac- 
tion. Earlier measurements of two-proton correlation 
functions in intermediate energy heavy ion collisions per- 
formed by Lynch et al. |9j provided the evidence for 
particle emission from localized highly-excited regions. 
After that, many experimental groups have investigated 
the momentum correlation functions in various aspects, 
such as for the unstable particle populations [T3 , the de- 
pendence on the impact parameter [TT| . the dependence 
of the total momentum of nucleon pairs |12| , the depen- 
dence of the isospin of the emitting source 0] and so 
on. More interestingly, HBT technique has been used to 
construct neutron-neutron correlation function which is 
useful to investigate the properties of neutron-halo nu- 
clei IH1 • The details about the nuclear Equation 
of State (EOS) and the collision dynamics could be re- 
vealed from the correlation function by the comparison 
between the experimental data and the transport model 
calculations. 

However, in most studies of p — p correlation functions, 
the HBT strength at 20 MeV/c of the proton-proton rel- 
ative momentum is taken a unique quantity to deter- 
mine the source size and / or emission time of two proton 
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emission. In a recent analysis for the HBT data below 
100 MeV/nucleon with the imaging method, Verde et 
al. show that the width of correlation function provides 
the information of the source size of the fast dynami- 
cal component while the peak of the correlation function 
is sensitive to relative yield from slow and fast emission 
components [rH Il8| . In addition, it was claimed that 
the proton emission from slow statistical component was 
not suitably treated in the conventional transport model, 
such as BUU model 19] . In light of the above studies, the 
whole shape of the correlation function is important for 
deducing space-time information of the emission source 
|17| . For heavy ion reactions with mid- heavy projectile 
and target combination in the Fermi energy domain, the 
slow emission component should not be neglected. To 
minimize the complication of the slow and fast compo- 
nent of proton emission in correlation function, it might 
be useful to choose light reaction systems at higher beam 
energies. In this context, we will use C + 12 C system to 
investigate momentum correlation functions above 100 
MeV/nucleon in this work. 

To understand the details of collisions for different re- 
actions by the HBT studies, a reliable simulation of the 
collision dynamical process for the heavy ion reaction is 
required. The simulation shall give a reasonable treat- 
ment of fragment formation after the final state inter- 
action. There are some good event-generator models 
to describe the collision process. In relativistic heavy 
ion collisions, successful models include both the string- 
hadronic-like models, such as relativistic quantum molec- 
ular dynamics and the parton cascade model etc. In in- 
termediate energy region, the successful transport model 
includes BUU model 12011 and Quantum Molecular Dy- 
namics (QMD) model [2l|. From those event-generator 
models one can obtain the phase space of the emitted 
particles with different parameters of the EOS and then 
construct the momentum correlation function. Recently, 
the nuclear symmetry energy dependence of HBT has 
been also explored through the isospin dependent BUU 
model [22^ in intermediate energy heavy ion collisions. 
Moreover, the nuclear binding energy and separation en- 
ergy dependences of the HBT strength have been inves- 
tigated with the help of the Isospin Dependent Quantum 
Molecular Dynamics model (IDQMD) [H 0]. In this 
paper, more features of momentum correlation function 
are reported by using the IDQMD model. 

The paper is organized as follows: in Sec. II we de- 
scribe the HBT technique and the IDQMD model. The 
stability of the IDQMD model is checked; Section III 
presents the results and discussions. We discuss the influ- 
ences of the following ingredients: initialization of projec- 
tile and target, emission time of nucleons, evolution time 
of the reaction, gate on the total momentum of nucleon- 
nucleon pair, soft and stiff momentum dependent EOS, 
in-medium nucleon-nucleon cross section, impact param- 
eter and incident energy etc. Finally the conclusions are 
presented in Section IV. 



II. HBT TECHNIQUE AND THE IDQMD 
MODEL 

Firstly, we recall the HBT technique. As we know, the 
wave function of relative motion of light identical parti- 
cles is modified by the final-state interaction and quan- 
tum statistical symmetries when they arc emitted in close 
proximity in space-time, and this is the principle of in- 
tensity interferometry, i.e. HBT. In standard Koonin- 
Pratt formalism |25t 1261 l27j , the two-particle correlation 
function is obtained by convoluting the emission function 
g(jp,x), i.e., the probability for emitting a particle with 
momentum p from the space-time point x = (r, t), with 
the relative wave function of two particles, 

r(T> _ fd 4 x 1 d 4 x 2 g(P/2,x 1 )g(P/2,x 2 )\cj ) ( q ,r)\ 2 
° l ' qj fd^x 1 g(P/2,x 1 )Jd^x 2 g(P/2,x 2 ) ' 

(1) 

where P(= pi + P2) and q(=g(pi — P2)) are the total 
and relative momenta of the particle pair respectively, 
and 0(q, r) is the relative two-particle wave function with 
their relative position r = (r 2 — ri)— |(vi + V2)(*2 — h)- 
This approach is very useful in studying effects of nuclear 
equation of state and nucleon-nucleon cross sections on 
the reaction dynamics of intermediate energy heavy-ion 
collisions p|. 

In the viewpoint of theoretical simulation, the cor- 
relation function can be established by using an event 
generator that produces the phase space information by 
modelling the collision dynamics and particle production. 
The event-generator correlation functions are then con- 
structed from the positions and momenta representing 
the single-particle emission distribution at the time of the 
last strong interaction, i.e. at freeze-out. In the present 
work, the event-generator is the Isospin-dependent Quan- 
tum Molecular Dynamics transport model |2l| . which 
has been successfully applied to the HBT studies of the 
heavy-ion collisions for neutron-rich nuclei induced reac- 
tions |23l 12 l| . Using a computation code named CRAB 
( Correlation After Burner ) of Pratt [23, which takes 
into account final-state nucleon-nucleon interactions, we 
have evaluated two-nucleon correlation functions from 
the emission function given by the IDQMD model. In 
the following part, we shall introduce the model briefly. 

The Quantum Molecular Dynamics approach is a many 
body theory to describe heavy ion reactions from inter- 
mediate energy to 2 GeV/nucleon (2^. It includes sev- 
eral important ingredients: initialization of the target 
and the projectile; nucleon propagation in the effective 
potential; nucleon-nucleon collisions in nuclear medium; 
Pauli blocking effect and the numerical test. A general 
review about the QMD model can be found in 21 1. The 
IDQMD model is based on the QMD model affiliating 
the isospin factors. 

The heavy-ion collision dynamics at intermediate ener- 
gies is mainly governed by three components: the mean 
field, two-body collisions and Pauli blocking. Therefore, 
for an isospin-dependent reaction dynamics model it is 
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important to affiliate isospin degrees of freedom with the 
above three components. In addition, the sampling of 
phase space of neutrons and protons in the initialization 
should be treated separately because of larger difference 
between neutron and proton density distributions for nu- 
clei far from the /3-stability line. For exotic neutron-rich 
nucleus one should sample a stable initialized nucleus 
with neutron-skin or halo structure so that one can di- 
rectly incorporate nuclear structure effects into a micro- 
scopic transport process. The IDQMD model has been 
improved based on the above ideas and the details will 
be given in the following. 

In the present calculations the interaction potential in 
the IDQMD is determined as follows: 

U(p) = [/Sky + yCoul +f/8 ym + ^Yuk +[/ MDI +[/ Pauli ; ^ 

where f/ Sky is the density-dependent Skyrme potential 
and it reads when the momentum dependent potential is 
included 

[/Sky = a( Z) + p(P_yt + t 4 ln 2 [e(^)2/3 + !j P. j (3) 
Po Po Po Po 

where p and p are total nucleon density and its normal 
value, respectively. The parameters a, (3, 7, t± and e 
are related to the nuclear equation of state and listed in 
Table. 1. V c is Coulomb potential. JJ Yuk is the Yukawa 
potential , 

U Yuk = t 3 eXP t l " ' (4) 
|ri— ra | 
m 

where m = 0.8 fm. [/ MDI is the momentum dependent 
interaction |29| . 

[/MDI = t4ln 2 [t5( ^_^| )2 + 1]j P. (5) 

PO 

where pi and p~2 are the momentum of two interacting 
nucleons. [/ Pauh is the Pauli potential , 

(6) 

where 

^ _ f 1 for neutron-neutron or proton-proton 
PiPi S for neutron-proton. 

The parameters V p , po and qo is 30 MeV, 400 MeV/c and 
5.64 fm, respectively. U sym is the symmetry potential. 
In the present calculation, we use U syirL = C sym ~ T a 
where C sym is the strength of symmetry potential, taking 
the value of 32 MeV. p n and p p are the neutron density 
and the proton density, respectively. r z is the zth com- 
ponent of the isospin degree of freedom, which equals 1 
or -f for neutrons or protons, respectively. The param- 
eters of the interaction potential are given in Table. 1 
where K — 200 or 380 MeV means the soft- or the stiff- 
momentum dependent potential, respectively. 
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Tabic 1. The parameters of the interaction potentials 

The in-medium nucleon-nucleon (NN) cross section can 
be parameterized as isospin dependent from the available 
experimental data. Studies of collective flow in HIC at 
intermediate energies revealed the reduction of the in- 
medium NN cross sections (3(1 El US- An empirical 
expression of the in-medium NN cross section [3l| is used: 

= (1 + /-)4n° (7) 
Po 

with the factor / w - 0.2 which has been found to better 
reproduce the flow data J30j . Here t^n ' s ^ ne experimen- 
tal NN cross section |33j . The free neutron-proton cross 
section cr^° is about a factor of 3 times larger than the 
free neutron-neutron or proton-proton cross section be- 
low about 400 MeV/nucleon in the laboratory energy. 
It should be mentioned that the relationship between 
the neutron-proton cross section and neutron-neutron 
(proton-proton) cross section depends also on the mod- 
ification of the nuclear density distributions during the 
reactions. 

The Pauli blocking effect in IDQMD model is treated 
separately for the neutron and the proton: Whenever a 
collision occurs, we assume that each nucleon occupies a 
six-dimensional sphere with a volume of h 3 /2 in the phase 
space (considering the spin degree of freedom) , and then 
calculate the phase volume, V, of the scattered nucleons 
being occupied by the rest nucleons with the same isospin 
as that of the scattered ones. We then compare 2V/h 3 
with a random number and decide whether the collision 
is blocked or not. 

When the initialization of the projectile and the target 
is taken in the IDQMD, the density distributions of pro- 
tons and neutrons are distinguished from each other. The 
references of neutron and proton density distributions 
for the initial projectile and target nuclei in IDQMD are 
taken from the Skyrme-Hartree-Fock (SHF) method with 
parameter set SKM*. Using this density distribution, we 
could get the initial coordinate of nucleons in nuclei in 
terms of the Monte-Carlo sampling method. The mo- 
mentum distribution of nucleons is generated by means 
of the local Fermi gas approximation: 

P* (f) = h(3ir 2 Pi (r))i, (i = n,p). (8) 
In the model, the radial density can be written as: 

x [exp(-±) - exp{—^)} (9) 
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where L is the so-called Gaussian wave width (here L = 
2.16 fm 2 ). 

Stability of the initialized nucleus has been checked by 
the time evolution of the system at zero temperature |24| . 
The accepted configurations are quite stable: only a few 
percentage nucleons escape from the nucleus till 200 fm/c 
in the intermediate energy domain. In addition, the sta- 
bility is also checked by tracking the time evolutions of 
the average binding energy and the root mean square ra- 
dius of the initialized nucleus and good enough stability 
has been found. Lighter nuclei are somewhat less sta- 
ble. One or two out of ten nuclei lose a nucleon in the 
required time span. To avoid taking unstable initializa- 
tion of projectile and target in the IDQMD calculation, 
we only select the initialization samples of nuclei which 
meet the required stability. 

Nuclear fragments are constructed by a modified 
isospin-dependent coalescence model, in which particles 
with relative momentum smaller than po = 300 MeV/c 
and relative distance smaller than Rq ~ 3.5 fm will be 
combined into a cluster. 

In our calculations, the reactions of C-isotopes with 
12 C have been performed. Most simulations are done in 
head-on collisions (6 = fm) at 100 MeV/nucleon or 800 
MeV/nucleon. The momentum correlation function are 
constructed by the phase space points at 200 fm/c when 
the system is basically at the freeze-out stage. 



III. RESULTS AND DISCUSSIONS 

A. Stable initialization vs random initialization 

Fig-fflshows proton- neutron correlation functions from 
the reactions induced by a chain of C-isotopes projectile 
on 12 C target at 800 MeV/nucleon of incident energy and 
head-on collisions (impact parameter 6 = fm) when the 
suitable selection of the stable initialization is used. In 
this figure the HBT strength of each isotope in lower 
relative momentum can be separated. If we plot this 
strength at 5 MeV/c (C p _„) as a function of the mean 
binding energy (Eb) of the projectiles (C-isotopes), we 
find that there exists an approximate linear relationship 
between C v - n and Eb as it is shown by the solid circles 
in Fig. El 

In order to investigate the importance of initialization 
in the HBT study, a comparison has been performed us- 
ing a random initialization for the projectile and the tar- 
get in the IDQMD based on Monte-Carlo sampling. In 
random initialization, the phase space of the nucleons is 
generated with Monte Carlo random sampling. During 
sampling of each nucleon, if the distance between every 
two nucleons is larger than 1.5 fm and the product of 
the space radius and the momentum radius between two 
nucleons is larger than one constant according to the un- 
certain relationship, the sampling will be accepted. In 
this case, there is no additional requirement of the bind- 
ing energy of the projectile or the target. 
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FIG. 1: (Color online) The proton-neutron correlation func- 
tion Cp-n for the reactions of a chain of C-isotopes with 12 C 
at 800 MeV/nucleon and b = fm using the stable initial- 
ization of projectile and target. The meaning of symbols are 
illustrated in the insert. 



Taken such a random initialization, the initial phase 
space of the nucleons for the projectile and the target 
is different event by event. Through the transport pro- 
cess of IDQMD, we can obtain the momentum correla- 
tion function and extract the strength of HBT in the final 
states. Fig. |3 shows the neutron-proton HBT for differ- 
ent C-isotopes by random sampled initialization. Similar 
to Fig. 01 the strength at 5 MeV/c is not a constant 
for different isotopes induced reactions which is shown 
by open circles in Fig. [21 However, two apparent dif- 
ferences can be observed between the results of stable 
initialization and random initialization. One is the dif- 
ference in the magnitude of the HBT strength. The val- 
ues of C p _„ at 5 MeV/c with the random initialization 
are less than those with the stable initialization. For the 
random initialization, the initial phase space which may 
not meet the requirement of the ground state as required 
in SHF calculation fluctuates event by event. In this case 
the tightness of initial nucleus becomes weaker than that 
with stable initialization of phase space which is sam- 
pled by the SHF density distribution. Thus, the values 
of Cp~ n at 5 MeV/c which can embody the tightness be- 
tween the nucleons become smaller. The other difference 
is the slope of C p - n vs Eb. A steeper linear relationship 
is observed for the stable initial phase space, while the 
dependence of C p - n vs Ej, becomes weaker in the random 
case. 

The comparison between both different initialization 
methods indicates that the reasonable initial phase space 
of the projectile which is sampled by the experimental Ef, 
and SHF density calculation in the IDQMD is important 
and suitable to investigate the dependence of the bind- 
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FIG. 2: (Color online) The relationship between the proton- 
neutron correlation function C p - n at 5 MeV/c and the aver- 
age binding energy per nucleon of C-isotopes. The filled cir- 
cles represent the results with the stable initial phase space 
taking in the IDQMD model according to the SHF density 
distribution. The open circles represent the results with the 
random initial phase space. The collisions were simulated at 
800 MeV/nucleon and b = fm. The target is 12 C. The lines 
are linear fits to guide the eyes. 



FIG. 3: (Color online) The proton-neutron correlation func- 
tion Cp-n for the reactions of a chain of C-isotopes with 12 C 
target at 800 MeV/nucleon using the random initialization of 
the projectile and the target. Symbols which represent the 
different isotopes are illustrated in the insert. 



ing energy for some observables. In the following calcula- 
tions, we will use the stable initialization to study some 
features of momentum correlation functions. 



B. Emission time of nucleons and evolution time 
of reaction 

In the intermediate energy domain, the extraction of 
the space-time information is further complicated by two 
effects. One is the presence of multiple sources of par- 
ticle emission [34]], another is the different time scale of 
statistical and dynamical emission from equilibrium and 
non-equilibrium sources [l7l lisl ] . The total momentum- 
gated correlation function can be used to investigate the 
later effects. Understanding the emission time sequence 
of neutrons and protons will be helpful to understand 
the nuclear interaction. It also might be sensitive to the 
nuclear equation of state. It has been shown that the 
emission times of the nucleons are related to their ki- 
netic energies. Generally, earlier emitted nucleons have 
higher energies than later ones. Some results have been 
demonstrated in the experiments |35j |. 

Fig. 01 shows the time evolution of the particle emis- 
sion rate by using the stiff momentum-dependent (Stiff 
M) EOS. The open circles connected with solid line show 
the proton emission rate and the filled ones represent the 
neutron emission rate. The collisions are simulated for 
18 C + 12 C at 100 MeV/nucleon and b = fm. In the fol- 
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FIG. 4: (Color online) The time evolution of the emission rate of 
nucleons for 18 C + C at 100 MeV/nucleon and 6 = fm with the 
different gates of the total nucleon-nucleon pair momentum. 
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lowing discussions, except the special note, the collisions 
are all performed in the above circumstance. 



From this figure, nucleons begin to be emitted around 
20 fm/c and their emission rates reach to their maximum 
15-20 fm/c later. With increasing total nucleon-nucleon 
momentum Ptot? peak of the emission rate of nucleons be- 
comes larger and its corresponding time tends to earlier 
time. This indicates that nucleons with higher total pair 
momentum are emitted earlier. Higher momentum nucle- 
ons mostly belong to pre-equilibrium emission nucleons 
and essentially originate from higher density regions. In 
contrast, lower momentum nucleons are mostly emitted 
from equilibrium-like sources. For neutron-proton pairs 
with lower total momentum, the emission rate of neutron 
and of proton is almost synchronous, ie. there is no obvi- 
ous difference of the emission sequence between neutron 
and proton. However, for the neutron-proton pair with 
larger total momentum, the emission rate of neutrons 
reaches to the peak value earlier than that of protons, it 
means that, on average, neutrons are emitted earlier than 
protons. The reason why the emission rate of neutrons is 
larger than that of protons stems from the neutron-rich 
content of projectile. Two interpretations seem to be 
possible. On one hand, the symmetry potential term in 
Eq. (2) plays an important role in controlling the emission 
of nucleons. In neutron-rich projectile induced reaction, 
protons could feel stronger attractive potential due to 
neighboring neutrons which results in more bound pro- 
tons for disassembling sources. Viceversa, neutrons will, 
on one side, feel the stronger repulsive interaction due to 
more neutron-neutron pairs and, on the other side, feel 
smaller attractive potential because of the decreasing of 
the average assortative number of nearest neighboring 
protons for a certain neutron for increasing isospin of the 
source. Both reasons will lead to produce more unbound 
neutrons for disassembling sources with higher isospin 

M. 



Experimentally, the momentum correlation function 
reflects the information in the final state of the reaction. 
Theoretically, the final state can be seen as the state at 
the freeze-out of system. Since the HBT is sensitive to 
the space-time information, we shall investigate the HBT 
at different evolution times of reaction. Fig. [S] shows the 
correlation functions of neutron-neutron, proton-neutron 
and proton-proton when the evolution time of reaction t 
= 150, 200, 250, 300 and 400 fm/c. Generally, HBT val- 
ues become smaller in later evolution time due to the 
weakness of the nucleon-nucleon correlation when the 
system is diluted. However, we can roughly say that 
the HBT values do not change dramatically after t — 
200 fm/c after t = 200 fm/c is compared to earlier times. 
Because of the limited computation resource, thereafter 
we investigate the features of HBT for all systems when 
t = 200 fm/c. 
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FIG. 5: (Color online) The momentum correlation functions of 
n-n, p-n and p-p pairs for 18 C + C at 100 MeV/nucleon and 6 = 
fm are constructed in different evolution time of the system. The 
evolution time is illustrated in the insert. 



C. Gate of the total momentum 

Since the magnitude of the total pair momentum is re- 
lated to the nucleon emission time, we shall discuss the 
effect of the total pair momentum on HBT in this section. 
Earlier emission time induces stronger correlation, larger 
total momentum thus contributes to the strength of cor- 
relation function too. We discuss the calculations with 
different total nucleon-nucleon pair momentum ( P to t ) 
in the following. The results are shown in Fig. In this 
figure, three types of nucleon-nucleon correlation func- 
tion, namely that of n-n, n-p and p-p, are shown. From 
the left panels, it is cleanly observed that the higher the 
Ptot, the higher the HBT value in lower relative momen- 
tum. While, from the right panels, it is observed that 
HBT value in the region of lower relative momentum is 
higher with the soft EOS than that with the stiff EOS. 

To see an overall trend of strength of HBT versus Ptot , 
we plot the value of C n _„ (C p _„) at 5 MeV/c or C p - P 
at 20 MeV/c as a function of P tot in Fig. □ The filled 
circles connected with the solid line present the results 
with the soft momentum-dependent EOS and the open 
circles connected with the solid line show the ones with 
the stiff momentum-dependent EOS. 

From the figure, it is clear that the strength of two- 
nucleon correlation function is smaller at lower total pair 
momentum than that of the higher one. As it is shown 
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FIG. 6: (Color online) The two-nucleon correlation function C n -n, 
C p -„ and C p - P for 18 C + C at 100 MeV/nucleon and b = fm with 
the different total pair momentum (Ptot) cuts or EOS parameters: 
the left panels correspond to the different gates of Ptot while the 
right panels correspond to different EOS for Ptot = 400 MeV/c. 



FIG. 8: (Color online) The comparison between the corre- 
lation function with soft momentum-dependent (Soft M) and 
stiff momentum-dependent (Stiff M) EOS for 18 C + C at 100 
MeV/nucleon and b = Ofm in IDQMD calculations. The left pan- 
els correspond to the cases of Ptot < 200 MeV/c while the right 
panels correspond to the cases of Ptot > 400 MeV/c. The filled 
and open circles represent the results with Soft- and Stiff- M EOS, 
respectively. 
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the nucleons with lower total momenta are emitted later 
than those with higher total momenta which naturally 
reduces the HBT strength. The calculated difference 
indicates the qualitative characterization of the emis- 
sion process during the collisions. On the other hand, 
the tendency is similar despite the different combina- 
tion of nucleon-nucleon pair. Experimental results for 
the momentum-gated nucleon pairs show a trend similar 
to the one that can be found in the literature, see for 
example Ref. |HD2. 

In addition, one can find that the HBT strength with 
the stiff momentum-dependent EOS is higher than that 
with the soft one. The influence of the different momen- 
tum dependent EOS on HBT will be discussed in the 
next section in details. 



D. Soft and stiff momentum-dependent potential 



FIG. 7: (Color online) The HBT strength of two-nucleon correla- 
tion function for n-n and p-n at 5 MeV/ c or at p-p at 20 MeV/ c for 
18 C + C at 100 MeV/nucleon as a function of the gate of the total 
momentum (Ptot) of the nucleon-nucleon pairs in the collisions. 



The EOS is considered an important property of nu- 
clear matter and several researches have been performed 
to investigate on the EOS of the finite nuclear matter 
[33. In this section, we will show the results of correla- 
tion functions with different EOS parameters. 

In previous studies, the role of different potentials, 
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i.e, the soft and stiff potential in the transport mod- 
els, has been investigated via some physical observables. 
However, in most case, the potential does not include 
the momentum dependent term, eg. in HBT studies 
with isospin dependent BUU j^. In this work we used 
a momentum-dependent part in the potential, namely 
the soft momentum-dependent potential (Soft M EOS) 
and the stiff momentum-dependent potential (Stiff M 
EOS), we use also an isospin-dependent potential and 
calculated nucleon-nucleon correlation functions with the 
above mentioned two types of momentum-dependent po- 
tentials. The results are shown in Fig. [S] 

In Fig. [3] the filled circles connected with solid line 
represent the correlation function with Soft M EOS and 
the open ones are the results with the one with Stiff M 
EOS. The left parts are the HBT results with P tot < 
200 MeV/c and the right parts are the results Ptot > 
400 MeV/ c. From the figure, it is clear that the correla- 
tion function with the stiff potential is larger than that 
with the soft potential, which is similar to the calculation 
results with BUU • Stiff potential makes the compres- 
sion of nucleonic matter difficult compared with the case 
with the soft one and leads to larger emission rate and 
earlier average emission time of nucleons, which leads to 
stronger correlation function for the stiff potential. 



E. In-medium nucleon-nucleon cross section 
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FIG. 9: (Color online) The two-nucleon correlation function Cnn> 
C p -„ and C p - P for 18 C + C at 100 MeV/nucleon and b = fm with 
different in-medium nucleon-nucleon cross section (^NN d /°NN ) 



The effect of the in-medium nucleon-nucleon cross sec- 
tion (o-™g d ) is discussed in this section. We use a value 
of 

°NN d which is different from the free nucleon-nucleon 
cross section (c^ e ) to investigate its influence on momen- 
tum correlation functions. Fig. shows NN, NP and PP 
correlation functions for different and EOS. Slightly 
larger values of the HBT is found for larger ct^n 1 ) espe- 
cially for 0.8 er<™°, as well as for the stiff EOS. Fig. EH 
shows the HBT strengths at 5 MeV/c (for NN and NP) 
or at 20 MeV/c (for PP) as a function of cr™ d in the 
stiff and soft EOS. The strength increases with <J^ d and 
with the stiffness of the EOS. This can be understood 
by the following arguments: with the increasing of the 
in-medium nucleon-nucleon cross section, the collision 
rate between two nucleons increases and consequently 
the system reaches to the equilibrium stage faster. Be- 
fore equilibrium is reached, more pre-equilibrium nucle- 
ons are emitted, which makes the strength of the corre- 
lation function larger. 



F. Impact parameter dependence 

Considering the importance of the Wigner function 
which depends on the impact parameter, we shall inves- 
tigate on momentum correlation function at different im- 
pact parameter. Collisions of 18 C + 12 C are performed 
at 100 MeV/nucleon and at the impact parameters of 0, 
1, 2, 3, 4 and 5 fm. The calculated total momentum inte- 




FIG. 10: (Color online) The HBT strengths of NN (a) and PN (b) 
at 5 MeV/ c and PP (c) at 20 MeV/ c as a function of the in-medium 
nucleon-nucleon cross section (^nnV^nn ) f rom the reaction 18 C 
+ 12 C at 100 MeV/nucleon and b = fm. The open and filled 
symbols represent stiff and soft momentum-dependent potential, 
respectively. 
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FIG. 12: (Color online) The strength of the nucleon-nucleon cor- 
relation functions of n-n (a) and NP (b) at 5 MeV/c or PP (C) 
at 20 MeV/c for 18 C + C at 100 MeV/nucleon as a function of 
impact parameter. The open and filled symbols represent stiff and 
soft momentum-dependent potential, respectively. 



FIG. 11: (Color online) Momentum correlation functions for n- 
n, p-n and p-p are calculated in different impact parameters (left 
panels) or with different EOS at fixed b = 3 fm for 18 C + C at 100 
MeV/nucleon. 



grated correlation functions are shown in Fig. 1111 Larger 
HBT values for central collisions or stiff EOS are pre- 
dicted. Their strengths are shown in Fig. 

In Ref. [llj , some explanations on the effect of impact 
parameter have been presented. As we already know, 
the strength of the correlation function mainly depends 
on the emission time and the source size. From the fig- 
ure, there exists a large difference between soft and stiff 
potentials. Secondly, the strength of correlation func- 
tion becomes weaker with increasing impact parameter 
on both soft and stiff momentum-dependent EOS. This 
indicates that the stiffness of the potential in the IDQMD 
does not change the tendency of the HBT strength of 
HBT with increasing impact parameter. On the other 
hand, the behavior of the HBT strength with impact pa- 
rameter might reflect the changed size of the emitting 
source. In central collisions, nucleon-nucleon collisions 
are very frequent and emitted nucleons are mostly from 
one compact and hot dense region. Therefore, the HBT 
strength is larger because of the smaller source size, if 
compared to peripheral collisions. Also some other rea- 
sons including Fermi jets may contribute to this effect 



G. Incident energy dependence 

The influence of the incident energy on HBT is inves- 
tigated in this section. Fig. 1131 shows the momentum 
correlation function at different beam energies (left pan- 
els) and different EOS (right panels). Fig.^Jdisplays the 
calculated HBT strengths of n-n (a) , p-n (b) and p-p (c) 
as a function of beam energy for head-on collisions. From 
the figure, both stiff and soft momentum potential pro- 
vides a similar evolution of the HBT strength with beam 
energy. Interestingly, the HBT strength first increases 
with the incident energy and reaches a peak around 100 
MeV/nucleon and then decreases at higher incident en- 
ergies. The raise of the HBT strength at lower incident 
energies can be essentially attributed to drastic nucleon- 
nucleon collision at higher energies, which results, on av- 
erage, in earlier emission of nucleons. However, the de- 
creasing behavior of HBT strengths with incident ener- 
gies higher than 100 McV/nuclcon cannot be explained 
in this way. 

In order to discuss the possible origin of the complex 
behavior of the HBT strength with incident energy, we 
investigate the maximum nucleon emission rate in the 
whole time evolution as a function of beam energy in 
Fig. 1151 From the figure, the maximum nucleon emis- 
sion rates show a rise and fall behavior similar to the 
one observed in Fig. ^| for the HBT strength. Around 
100 MeV/nucleon, the nucleon emission rate is maxi- 
mum. This energy coincides with the peak position of 
the HBT strength in Fig. From this point of view, 
the HBT strength could be correlated to the maximum 
nucleon emission rate during heavy ion collisions. The 
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FIG. 13: (Color online) The momentum correlation function of n-n 
(a), p-n (b) and p-p (c) pairs in different beam energies in head-on 
collisions for soft M EOS (left panels) or different EOS at Ebeam 
= 200 MeV/nucleon (right panels). 



higher the nucleon emission rate, the stronger the HBT 
strength. Of course, another possible explanation for the 
falling branch of HBT strength can be attributed to the 
source size in the final state (t = 200 fm/c). It is ex- 
pected that in the energy range of a few hundreds of 
MeV/nucleon, where the repulsive nucleon-nucleon inter- 
action is dominant, a more dilute system might develop 
after freeze-out. So the fall of HBT with beam energy 
above 100 MeV/nucleon indicates that the source size 
has been much expanded in the final state at higher en- 
ergies. In contrast, the attractive mean field competes 
with the repulsive nucleon-nucleon interaction below 100 
MeV/nucleon. Shorter emission times with increasing 
beam energies may play a dominant role in determining 
the behavior of the HBT strength. Fig.^|also illustrates 
illustrates that the emission rate of neutrons is larger 
than that of protons and that the neutrons are emitted 
earlier than protons. 



IV. SUMMARY 

In summary, isospin-dependent quantum dynamics 
model has been used as an event-generator to study 
the momentum correlation function of neutron-neutron, 
proton-proton and neutron-proton correlation functions 
for C + 12 C. An approximate linear relationship between 
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the strength of correlation function at low relative mo- 
mentum and the mean binding energy of projectiles has 
been revealed and the influence of the different initial- 
ization methods of the projectiles on momentum corre- 
lation function has been also investigated. In addition, 
some different physical ingredients which have effects on 
momentum correlation functions have been shown in this 
work. These ingredients include nuclear equation of state 
and in-medium nucleon-nucleon cross section. It shows 
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that the stiff momentum-dependent EOS or a larger in- 
medium nucleon-nucleon cross section results in stronger 
HBT strength than the soft momentum-dependent EOS 
or a smaller in-medium nucleon-nucleon cross section. 
Some other aspects which have also effects on momen- 
tum correlation function, including the time evolution of 
the reaction system, the impact parameter, the gate of 
the total momenta of nucleon-nucleon pairs and the in- 
cident energy are also explored. Results show that the 
HBT strength decreases with increasing impact parame- 
ter due to the earlier emission or a compact source size 
in central collisions. A positive correlation between the 
gate of total nucleon-nucleon pair momentum and the 
HBT strength can be explained by the earlier emission 
time of nucleons with higher total momentum. The rise 
and fall of the HBT strengths with the beam energy can 
be interpreted either in terms of the maximum nucleon 



emission rates during heavy-ion collisions or in terms of 
the dominant shorter emission time with increasing beam 
energy below 100 MeV/nucleon and the larger later-stage 
expansion of the source size with the increasing of beam 
energy above 100 MeV/nucleon. 
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